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INTRODUCTION 

Edge-flames play an important role in a number of microgravity investigations, and in the general 
study of flames. Examples include the candle-flame experiments carried out on board both the 
Space-Shuttle and the Mir Space Station, e.g.[l]; the flame-spread-over-liquid work carried out by 
H.Ross and W.Sirignano, amongst others, e.g. [2,3]; and lifted turbulent diffusion flames. In all 
of these configurations a local tw o-dimensionalflame structure can be identified which looks like 
a flame-slieet with an edge, and these structures exhibit dynamical behavior which characterizes 
them and distinguishes them from ad ho c2D flame structures. 

Edge-flames can exist in both a non-premixed context (edges of diffusion flames) and in a 
premixed context (edges of deflagrations), but the w ork reported here deals with the edges of 
diffusion flames. It is particularly relevant, we believe, to oscillations that ha vdseen seen in both 
the candle-flame context, and the flame-spread-over-liquid cortext. These oscillations are periodic 
edge-oscillations (in an appropriate reference frame), sans oscillation of the trailing diffusion flame. 

References [4,5] examine a simple model of an anchored edge-flame, and construct solutions 
using DNS. It is shown that if the Lewis number of the fuel is sufficiently large (the Lewis number 
of the oxidizer is taken to be 1), and the Damkohler number is sufficiently small, oscillating-edge 
solutions can be found. Oscillations are encouraged by an on-edge convectixe flow and the insertion 
of a cold probe, discouraged by an off-edge convectiv tflo w. 

In the present work, the nature of these oscillations is examined in more depth, using a variet y 
of numerical strategies. 
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Figure 1: (i) The model configuration (the low' eiand upper boundary's are at y = 0, y - 1); (ii) 
Response diagram for the steady 2D solutions. 
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Figure 2: (i) Unstable solution for D = 1.7 x 10 7 ; (ii)V ariationsof the period of oscillation with 
Damkohler number. 


THE CONFIGURATION AND STEADY SOLUTIONS. 


The model configuration, first introduced in [4], is sketched in Figure l(i). The discontinuity in 
data at the fuel-supply boundary introduces a weak anchor for an edge-flame nominally located in 
the region x > 0. The field equations are 

| (T, X,Y) = v 2 (t, ±-X, jLy) + (,, -or) DXYe~‘l T , (1) 

and these describe a ID frozen solution at x — oc, a ID diffusion flame solution at x — > -Hoc. 
The latter is characterized by the familiar S-shaped response of diffusion flames, with an upper 
strong-burning branch when D > D m i n . When 


11 2 

Q A- - 2 ’ = 2’ Lex ~ Ley - 2 > 9 = 1 , 0-6, T w = — 

(values that we shall use throughout), we find D m i n = 0.79 x 10 7 . 

Solutions of the 2D problem defined by (1) are constructed using 4th-order spatial differencing 
and Runge-Kutta time integration. F or some values of D, stationary solutions are unstable. T o 
calculate these using the unsteady code we fix the integral 


In 



DXYe~ 0/T . 


( 2 ) 


Then at eac h fractional time step in the R-K in tegration in which A", Y, and T are adjusted, D 
is also adjusted to fix /$>• The solutions constructed in this way aE c haracterized by the response 
of Figure l(ii) where the total heat output (ordinate) is Iq, since here q = 1. This reveals a 
detachment Damkohler number, D^t = 1.44 x 10 7 , below which 2D solutions do not exist. 

Unsteady DNS at. fixed D identifies stable steady solutions at late time if D > D ns , and 
oscillating solutions for some range of D in D < D ns . Figure 2(i) shows temperature variations 
at a fixed point when D = 1.7 x 10 7 , a valueless than D ns ; indeed D ns lies somewhere betw een 
1.9273 x 10 7 and 1.9580 x 10 7 . 

At this point we have identified three significant Damkohler numbers: the ID quenching value 
Dmin — 0.79 x 10 7 ; the 2D stationary- edge detachment value Ddct = 1-44 x 10 7 ; and the 2D 
neutral stability value D ns k, 1.94 x 10 7 . A fourth significant nunber is the minimum value for 
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which oscillating (non-detached, non-quenched) solutions can be obtained. This is approximately 
D = 1.583 x 10 7 , a value greater than D det . As D\ 1.583 x 10 7 , the amplitude of the oscillations 
increases, so that the maximum distance between the edge and the anchor point at (0, 0) increases. 
Apparently botlthe amplitude and the period of the oscillations approach infinity , albeit slowly. 
Thus Figure 2 (ii) shows the period for various D. The minimum D here is 1.5835 x 10 7 and x max 
for the computational domain is 8, to prevent failure due to disappearance of the flame from the 
domain. 


LINEAR STABILITY 


We have seen how steady 2D edge-flame solutions can be constructed for all D greater than a 
critical detachment value, whether they are stable or unstable. These solutions can be used as the 
foundation of a linear stability analysis. 

If w einclude 2 derivativ esin the description (1) and linearize about the steady-state solution 
by writing 

T = T s + T', etc. , (3) 


the perturbation variables satisfy the linear system 


d_ 

dt 


(X',Y',T') 


(£*£-•■) 




+ (-a x ,- ay ,q)D[xYe- d ' T ] , 

( 4 ) 


where we have assumed that the perturbations are proportional toe lkz , corresponding to corruga- 
tions in the z direction of specified wave-n uiher. The boundary conditions are 


T' = X' = Y' = 0 at 3/ = 0,1; 


d 

— ( )' -» 0 as |a:| -» oc. 
ox 


( 5 ) 


For the parameter values examined here, the relevant ID solution is stable for all D greater than 
the ID quenching value. 

A common strategy for discussing the solutions to systems sue has (4), (5) is to make the 
replacement ^ A and solv etlie eigenvalue problem by finite differencing and matrix methods. 
But an alternative is to use the unsteady solv erthat w eha vaised for all of the other numerical 
components of the problem. Time integration of (4), (5) with arbitrarily chosen initial conditions 
leads to an ever-grewing solution if D lies to the left of D ns . This growth can be monitored for 
one of the field variables, at a single point. Provided there is only a single unstable mode, both the 
gro wth rat.c7vf( A) and the frequency Im(A) can be deduced from this output. At the same time, 
a late-time snapshot of the field variables identifies the eigenfunction, e.g. Figure 3(i). The dipole 
nature of the reaction-rate eigenfunction, Figure 3(i), is simply a reflection of the oscillating nature 
of the disturbance. When the reaction center, the core of the edge, moves to the left, the reaction 
rate is augmented at points on the left, diminished at points on the right, as here. Note that the 
perturbation is confined to the edge, and does not affect the trailing diffusion flame. 

Figure 3(ii) shows variations in the growth rate Ke( A) with waw-number k when D = 1.777 x 
10 7 . V ariationsof frequency with k are modest, and w edo not sho wtliem. These results are 
typical of those that w eha vefound when D is in the unstable range. Further details and 3D 
DNS calculations that we have carried out will be reported elsewhere, but we note that the growth 
rate is a maximum when k — 0 and, consistent with this, the 3D DNS calculations show that 3D 
disturbances are suppressed in favor of 2D disturbances. 


NASA/CP— 2001-210826 


343 





Figure 3: (i) P erturbation reaction-rate colours: the left, group are positive, values 10 50 150 300 
450 600 750; the righ tgroup are negative, values -10-50 -150 -250 -300; (ii) V ariationsof growth 
rate with wa\e-number, D = 1.777 x 10 7 . 

CONCLUDING REMARKS 

In this study we have gained further insight into the nature of oscillating edge-flame solutions. But 
of greater interest, perhaps, is the demonstration that unsteady DNS can be used to construct 
unstable steady solutions and to solv ethe corresponding linear stability problem. The n umerical 
examination of stability problems is commonplace in the fluid mechanics literature, but rare in the 
combustion literature, and perhaps the results presented here will encourage others to pursue this 
strategy in a variet y of applications. 
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